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Abstract. In former papers we showed that during the decay of a neutron star's magnetic field under the influence 
of the Hall-drift, an unstable rise of small-scale field structures at the expense of the large-scale background field 
may happen. This linear stability analysis was based on the assumption of a uniform density throughout the 
neutron star crust, whereas in reality the density and all transport coefficients vary by many orders of magnitude. 
Here, we extend the investigation of the Hall-drift induced instability by considering realistic profiles of density 
and chemical composition, as well as background fields with more justified radial profiles. Two neutron star models 
are considered differing primarily in the assumption on the core matter equation of state. For their cooling history 
and radial profiles of density and composition we use known results to infer the conductivity profiles. These 
were fed into linear calculations of a dipolar field decay starting from various initial configurations. At different 
stages of the decay, snapshots of the magnetic fields at the equator were taken to yield background field profiles 
for the stability analysis. The main result is that the Hall instability may really occur in neutron star crusts. 
Characteristic growth times are in the order of < 10 4 . . . 10 6 yrs depending on cooling age and background field 
strength. The influence of the equation of state and of the initial field configuration is discussed. 
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1. Introduction 

Neutron stars (NSs) are carriers of the strongest magnetic 
fields which occur in nature. But, astonishingly enough, 
most of the quantitative studies of magnetic field decay 
in NS crusts consider only the linear induction equa- 
tion, i.e., a fi eld decay caused solely by Ohmic dissipa- 
tion (see, e.g..lUrpm et al.lll994UUrpin fc Konenkovll997t 
IPaee et al.ll2000j) . 

lYakovlev fc Shalvbkovl <ll99lli showed that in a two- 
component plasma the resistivity components parallel and 
perpendicular to the magnetic field coincide, as a result 
of which in turn the ambipolar drift disappears. This re- 
sult is applicable to the electrons in t he fully ionized crys- 
tal m atter of the NS crust, too (see lUrpin fc Shalvbkovl 
Il999lh Further, convcctivc motions will not exist since 
the crust is almost completely crystallized after, say, 10 4 
yrs. Therefore, the magnetic field evolution in the crust 
is solely determined by Ohmic diffusion and the so-called 
Hall-drift where the latter is the only non-linearity in 
this process (if the weak and therefore weakly nonlinear 
back-reaction of the magnetic field on the conductivity 
ten sor via Joule heating is discounted). Recently, two of 
us I Rheinhardt fc Geppertl 120021 iGeppert fc Rheinhardt! 



showed that for a (large-scale) magnetic back- 
ground field characterized by a sufficiently curved radial 
profile beyond a certain marginal field strength, the Hall- 
drift may cause an unstable growth of small-scale per- 
turbations. For a homogeneous medium these are concen- 
trated towards the boundary adjacent to vacuum. 

We performed this stability analysis based on the lin- 
earized induction equation with Hall-drift for a homoge- 
neous plane slab of finite thickness, infinitely extended in 
both horizontal dimensions and bound by vacuum and a 
medium of infinite conductivity at its upper and lower 
sides, respectively. While this model is perhaps an accept- 
able approximation of the NS crust with respect to geom- 
etry, the assumption of a spatially uniform density and 
chemical composition, which result in a uniform scalar 
conductivity and Hall-drift coefficient, is surely a very 
crude one. Actually, the density and thus the coefficients 
depending on it may vary thr oughout the crus t by many 
orders of magnitude (see, e.g.. |Paee et al-lEoOO]) . Less pro- 
nouncedly, the chemical composition (that is, mass num- 
ber A and atomic number Z) varies, too. Additionally, 
the scalar conductivity and the Hall-drift coefficient are 
dependent on the temperature and the impurity concen- 
tration. 
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The effect of the Hall-drift on the magnetic field 
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which includes quantitatively the effect of the crustal 
density stratification in this context is the one by 
IVainshtein et alJ l(200fll) . They showed that due to a den- 
sity gradient, the Hall-drift may create current sheets, 
where the field can be dissipated very efficiently. However, 
the occurrence of the mentioned Hall-instability could 
not be detected by them because in their linear analy- 
sis they considered a uniform background field only, i.e., 
a field which will not become unstable regardless which 
strength it has. In considering the non linear evolution 
of the toroidal field they neglected the coupling with the 
poloidal one, thus making an instability impossible either, 
see (0 below. 

Observational evidence for a decay of the large- 
scale (dipolar) magnetic field of middle aged pul- 
sars, being drastically accelerated in comparison with 
the purely Ohmic d ecay, has been discussed in 
IGeppert fe Rheinhardtl l)2002|) . Our reasoning there is 
based o n the detection of braking indices greater than 
three bv l.Tohnston fc Gallowavl ( 19991) , who used measure- 
ments of the rotational period P, and its temporal deriva- 
tive P, dating from different observational epochs. As a 
typical result we found that the decay times inferred from 
P and P may be smaller than 10 -4 times the (estimated) 
Ohmic decay time. A possible explanation for such a rapid 
field decay may be the Hall-instability which drains mag- 
netic energy out of the dipolar field and uses it for the 
build up of small-scale magnetic field structures, which 
eventually decay by Ohmic dissipation. 

The occurrence of small-scale magnetic field struc- 
tures in the crustal layers of NSs may cause further effects 
which are potentially observable. Beyond the decrease of 
the braking power of the magnetic field, the increased 
small-scale Lorentz forces can trigger a cracking of the 
crystallized crust. Also, due to a rapid and spatially con- 
centrated Ohmic decay of the small-scale field structures, 
hot spots may appear in t he surface layers. R ecent obser- 
vations, bot h in the X-ray jBecker et al.l2003|) . and in the 
radio range l|Gil et al.ll2002|) . support the idea of the ex- 
istence of strong small-scale magnetic field structures at 
the NS surface. 

Of course, none of the mentioned phenomena can 
be unambiguously attributed to the effect of the Hall- 
instabilit y. Some possible alternative exp lanations are dis- 
cussed in iGerroert fc Rheinhardtl (|2002|) . 

However, the simultaneous existence of a sufficient cur- 
vature of the background field profile and of a sufficiently 
large magnetization parameter cjsr e related to that field 
(where lobt c > 1 is a necessary condition) will very likely 
lead to the appearance of the Hall-instability. How vig- 
orously it develops for a realistically modelled crust and 
at which strength it saturates, remain as open questions. 



Also unanswered is, how far it may be hidden by other ef- 
fects of the field evolution as, e.g., the "normal" Hall cas- 
cade. The first of these questions will be addressed here. 

One of the most convincing indicators for the impor- 
tance of the Hall-drift in the crust is the evolution of 
the magnetization par a meter u>BT e there. In Fig. 1 of 
IGeppert fc Rheinhardtl l|2002|) this quantity is shown for 
different temperatures, that is, for different ages as a func- 
tion of density. These results are based on linear field de- 
cay calculations for a standard NS mod el with a medium 
EOS as presented bv lPaeeet all <|2000() . Clearly, as soon 
as the magnetic field strength exceeds 10 12 G, ujbt > 1 
in some regions of the crust, and the Hall-drift begins to 
dominate the Ohmic decay. 

In order to answer the question whether at all and 
how intensive the Hall-instability occurs in real NS crusts, 
we will consider crustal density profiles resulting from 
NS models based on stiff and medium equations of state 
(EOSs) of the core matter. As background magnetic fields 
which must - in comparison with the unstable field pertur- 
bations - evolve very sl owly, we use those , calculated by 
the methods described in lPaee et al.ll|200(t . The evolution 
of those large-scale fields is affected only by the density 
profile of the NS crust and its cooling history, both deter- 
mined essentially by the EOS. Furthermore, the crustal 
field evolution is affected by the chemical composition 
and impurity concentration within the crust, as well as 
by the initial strength and structure of the field which, in 
turn, reflect the processes at birth of the NS. These results 
are insofar incomplete, as they just do not take into ac- 
count the very efficient drain of magnetic energy from the 
background field during the period of the Hall-instability. 
They should, however, provide hints under which condi- 
tions, more realistic than those considered up to now, an 
episode of strongly non-linear magnetic field decay may 
take place. 

This paper follows the line s of thought as presented by 
iRheinhardt fc Geppert! l|2002() . and is organized as follows: 
In the next section the basic equations of the model are 
derived, and the assumed properties of the crustal matter 
are described together with the background field profiles 
we used. Section briefly outlines the method of solving 
the eigenvalue problem and provides and discusses the nu- 
merical results. Conclusions are drawn in Sect. 0] with a 
focus on possible observational consequences of the Hall- 
instability. 

2. Description of the model 

2.1. Basic equations, geometry and boundary 
conditions 

In the absence of convective motions and of ambipolar 
diffusion, conditions found in the solid crust of NSs, the 
equations governing the magnetic field evolution read 

B = — c curl | — - — (curlB + ujBT e (cm\B x e#)) ) 

V47TCT V V (1) 

div B = , 
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where r e is the relaxation time of the electrons and ll>b = 
e\B\/m* e c the electron Larmor frequency, with e being 
the elementary charge, m* the effective mass of an elec- 
tron, and c the speed of light, eg is the unit vector in 
B-direction. The scalar conductivity ctq is given by 



Ze 2 p{l-X n ) 



A 



m u rrr e 



(2) 



Here p denotes the (depth depending) density of the 
crustal matter, and m u the atomic mass unit. By virtue 
of their direct dependences on density, both m* and r e 
depend on the crustal depth, too. The mass and atomic 
numbers A and Z, respectively, also depend on depth 
as well as the fraction of free neutrons X„ (see , e.g., 
iNegele fc Vautherinl Il973l lHaensel fc Zdunikl llQQOft . Eq. 
Q lets one expect that the Hall-drift described by the 
term wsr e (curlB x ej) may become important for the 
magnetic field decay only if tosTe > 1. 

Linearization of Q with respect to a background field 
(= reference state) Bo yields: 

b = — curl (ji curl 6) 

— curl {a ( curl B x b + curl b x B )) (3) 
div b = 

describing the behavior of small perturbations b of the 
reference state. Here, the magnetic diffusivity rj, and the 
Hall coefficient a are given by 



n 



A 



47T(7n 



n- 



\B\ 



ZAnep{l-X n ) 



(4) 



Note, that the Hall parameter a is time-independent for a 
non-accreting NS while rj depends strongly on the crustal 
temperature via r e and, therefore, on time. In general, 
Bq is a decaying field, i.e., time-dependent, too. For the 
instability analysis we nevertheless want to treat them as 
constants. Then the following necessary condition has to 
be satisfied to take the results for valid: The growth time 
of an unstable mode has to be significantly shorter than 
all the characteristic times of changes of coefficients in 
iJSJ especially significantly shorter than the decay time of 
the background field. In the following this restriction will 
be referred to as the "background dynamics permissibility 
condition" . 

Employing the argume n ts from dynamo theory quoted 
in iRheinhardt fc Oeppertl l|2002l) , one can conclude that 
a curl Bq must not be a homogeneous field. More precisely, 
when interpreting this term as a velocity field, it must not 
be a rigid body motion, in order to be capable of delivering 
energy to b. Hence, in contrast to the homogeneous density 
model, now any inhomogeneous Bq may potentially enable 
the instability. 

Let us now specify the geometry of our model. 
Idealizing the spherical shell of the NS crust, we consider 
a plane slab which is infinitely extended both into the 
x- and y-directions, but has a finite thickness d in z— 
direction. We identify z with the crustal depth being zero 



at the surface of the NS. The background field is assumed 
to be parallel to the surfaces of the slab pointing, say, in 
a;-direction and to be dependent on the depth z only, i.e., 



B = B e x = f(z)e x 



(5) 



where in comparison to the homogeneous density model 
the minimum demand on f(z) can now be relaxed from 
quadratic to linear in z, if only a is not a constant (see 
above). Note, that due to this assumption a curl Bo x Bq 
represents a gradient. Thus the unperturbed evolution of 
the background field is not at all affected by the Hall-drift; 
in the absence of any electromotive force it would decay 
purely ohmically. 

XXX The choice (jSJ is not only motivated by the sake 
of simplicity, but to some extent justified by the physical 
conditions during the proto-NS stage (after the end of 
a possible field-generating phase): As long as the layers 
which later on form the crust are liquid, the magnetic 
field adjusts to be close to a magnetostatic equilibrium. 
Then it approximately obeys the condition 



curl I — curl B x B ) = 



with 



(6) 



being the elec tron number density, (see 



iThompson fe Duncan! Il993l Sect. 14.1) 1 . As the Hall- 
coef ficient a can be written in the fo rm a = c/Aim e e (see, 
e.g.. lOoldreich fc Reiseneggerlll992l) this is at the same 
time the condition for the Hall e.m.f. in Q to be a gra- 
dient and thus ineffective. However, as the modes of free 
decay in the crystallized crust are in general violating © , 
the crustal field will tend to deviate increasingly from the 
magnetostatic configuration. But at least for early stages, 
one may suppose that although ubt & might already be 
bigger than unity, the effect of the Hall term on the back- 
ground field is dominated by ohmic dissipation. XXX 

In Fig. Q the model geometry is shown with three 
different depth profiles of Bq, depicting qualitatively some 
of those we employed. They will be specified later. 

We decompose a perturbation b into poloidal and 
toroidal constituents, b = b p + b t , which can be repre- 
sented by scalar functions S and T, respectively, by virtue 
of the definitions 



b„ 



curl ( 



xVS). 



-e 2 xVT, (7) 



ensuring div 6 = for arbitrary S, T. 

For the sake of simplicity, we will confine ourselves to 
the study of plane wave solutions with respect to the x— 
and y-directions, thus making the ansatz 



(x,t) 



(z) exp (ikx + pr) , 



(8) 



where r denotes the time, k = (k x , k y ), x = (x,y) and p 
is a complex time increment. 

1 There exist even dipolar fields in spherical shells satisfying 
© with n e = n e (r) exactly. 



4 



Rheinhardt, Konenkov & Geppert: Hall-Instability in Neutron stars 




z = d 



Fig. 1. Sketch of the model geometry (Bq - background 
fields). In z > d a perfect conductor, in z < vacuum is 
assumed. 

Inserting @ with JHJ into J3J), we finally obtain two 
coupled ordinary differential equations for the scalars s 
and t: 

pt - r)(t" - k 2 t) - rj't' +iak x (f"s + k 2 fs - fs") 

+ ia'(k y ft + k x f's) = (9) 
ps - r](s" - k 2 s) +ia(k x ft - kyf's) = , 

where the dash denotes the derivative with respect to z, 
and k 2 = k 2 + k 2 . In comparison with the correspond- 
ing Eqs. (6) in iRheinhardt fc Geppertl l|2002|h the terms 
oc if and a' occur additionally, but they don't affect the 
energetic conclusions relevant for the existence of the in- 
stability. 

When completed with appropriate boundary condi- 
tions, these equations define an eigenvalue problem with 
respect to the time increment p. The boundary conditions 
chosen here are transition to vacuum at z = 0, and to a 
perfect conductor at z = d, respectively, mimicking the 
superconducting core below the bottom of the crust and 
an atmosphere with low conductivity above its surface. 
They read 



ks = t 



= 



s = rjt — ik y af t = 



for 
for 



2 = 

z = d 



(10) 



These conditions are equivalent to the requirements that 
all components of the magnetic field arc continuous across 
the vacuum boundary, and that neither the normal com- 
ponent of the magnetic field nor the tangential compo- 
ne nts of the electric field penetra te the core. For details 
see lGeppert fc Rheinhardtl l|2002l) . 

The signs of the wavenumbers are irrelevant for the 
eigenvalues of I© with (fTTTf) . since the transformations 



k x ► k x . 



p, s 
and 



ky > ky , P ► P 



t -> ~t 
(11) 



hold (* is the complex conjugate). Thus it is sufficient 
to consider the quadrant k Xl k v > of the fc^-fc^-plane 
only. On the basis of (|1 1|) , it can be easily concluded that 
changing the sign of / is irrelevant, too, because 



/-►-/, P 



t -> t* 



(12) 



As can be inferred from @ using standard energy 
arguments all solutions for k x = are damped oscil- 
lations, being either purely toroidal or purely poloidal. 



However, for k y = growing solutions with their poloidal 
and toroidal parts mutually coupled seem well possible if 
%{k x J*(af/ V ) (s*t)dz}>0. 



2.2. Crustal matter properties 

The main goal of the present paper is to enhance our 
knowledge on the conditions under which the Hall- 
instability occurs. As one step in this direction we con- 
sider a stratified crust characterized by a realistic density 
profile instead of a homogeneous crust. 





FP 


PS 


EOS 


medium 


stiff 


compactness GM/c 2 R 


0.2 


0.13 


central density (gem -3 ) 


1.17 ■ 10 15 


3.64 ■ 10 14 


NS radius (km) 


10.64 


16.4 


crust thickness (m) 


730 


3860 




Friedman & 


Pandharipande 


reference 


Pandharipande, 


& Smith, 




1981 


1975 



Table 1. Properties of the Pandharipande Smith (PS) 
and Friedman-Pandharipandc (FP) NS models; Mjvs = 
1.4M . 

The resistivity tensor is strongly dependent on den- 
sity, temperature, and chemical composition of the 
crust (see Eqs. (01, (@J; r e , as far as collisions with 
pho nons are concerned dep ends on temperature, too, 
see lUrpin fc Yakovlevl l|l980|) '). In turn, these quantities 
are in their temporal behavior mainly determined by 
the EOS of the NS matter, especially that of the core. 
Therefore, one has first to make a choic e among the mul - 
titude of proposed core EOSs (see, e.g.. Ivan RroerlllP.9lh . 
We decided to consider two of them, which are gener- 
ally accepted to be typical repres entatives of a medium 
jFriedman fc Pandharipanddll98l[ henceforth FP) and a 
stiff (|Pandharipande fc Smithll975l henceforth PS) EOS, 
respectively, thus probably covering an essential part of 
the obse rved NSs. (The ext reme soft EOSs, e.g., that de- 
rived bv lBavm et all l|l97l() are now conside red unlikely 
to be realized in nature, see IPage et alJkoOOj) Further on, 
we specify the crust in the FP case to consist of cold 
catalyzed matter, and in the PS case to be dominated 
by matter accreted and processed in the past. The cor- 
responding EOSs of the crus tal matter and its chemica l 
com positions were derived in iNegele fc Vautherinl l|l973l) 
and lHaensel fc Zdunikl ljl990l) . respectively. Since these 
two models are referred to and utilized in numerous pa- 
pers, we thus hope to facilitate discussion and application 
of our results. 

Some characteristic properties of the selected mod- 
els resulting from their density structure specified for 
Mns = 1.4M Q , are summarized in Table ^ Figure [2 
shows the density profiles. Additionally, the profiles of the 
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2x10 4x10 
depth (cm) 




depth (cm) 



Fig. 2. Radial density (solid) profiles for the FP (top) and 
PS (bottom) models. Additionally, for the latter the chem- 
ical composition, A, Z, (dashed and dash-dotted, respec- 
tively) is given. 



chemical composition, that is, A and Z, for the PS m odel 
as derived from the data in lHaensel fc Zdunikl l)l99fi) arc 
given. There exist a number of discrete depths (i.e., den- 
sities and pressures) where the preferred species of nuclei 
change almost abruptly causing discontinuities in A and 
Z. For the FP model, the corresponding data were given 
only implicitly in the form of a smooth fitting function for 
p{z). In both cases the crust-core boundary was assumed 
to be at a density of 2 x 10 14 gcm~ 3 whereas the crust 
surface was defined by the density 10 10 gcm~ 3 . 

The most important of the EOS's direct consequences 
is the degree of compactness: The smaller compactness 
(i.e., larger radius and crust thickness) of the PS model 
compared with the FP model of equal mass, results from 
its stiffcr EOS. Correspondingly, the latter model cools 
down significantly slower than the former for ages higher 
than 10 5 yrs, since a large compact ness (strong gravity ) 
inhibits the heat transfer (see, e.g.. iMisner et aflll973() . 
Fig. |3 sho ws the cool i ng cu rves for both models as cal- 
culated by Ivan Riper! (Il99lh . Throughout this paper we 
confine ourselves to NSs old enough (say, some 10 3 yrs) 
that the crustal temperature can be considered uniform, 
and the density profile is no longer changing in time. 
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- 
















J) 






\ 





10 10 
age (yr) 



Fig. 3. Cooling curves for the FP (squares) and PS (tri- 
angles) models. The symbols denote the ages at which the 
stability analysis is performed. 



Finally, Fig. Ogives the profiles of the magnetic diffu- 
sivity r\ and the Hall parameter a for both models. For 77 
different profiles are shown corresponding to the ages at 
which the stability analysis will be performed, but only 
one profile for a since it doesn't depend on temperature 
and hence not on age. For rj an impurity concentration of 
1% was assumed. Note, that in contrast to the PS model 
merely the implicit influence of A and Z on rj and a via the 
density is included for the FP model, the explicit one being 
neglected as well as the influence of X n . Since the cooling 
of a NS proceeds faster in the PS model, further cooling 
will not affect the conductivity after about 5 x 10 5 yrs in 
this case. Then phonons are practically no longer excited 
in the crustal crystal and the conductivity is determined 
by electron- impurity collisions alone. Hence, the diffusiv- 
ity doesn't change any more, while for the FP-model a 
significant number of phonons may be excited up to an 
age of approximately 5 x 10 6 yrs; therefore r\ depends on 
time until this age. 

According to the discontinuous behavior of A and Z in 
the PS model (see Fig. 0), a non-smooth behavior of its 
conductivity parameters appears in Fig.0] We decided not 
to smooth out these discontinuities because it had meant 
introducing unnecessarily some further arbitrariness into 
the model. 

2.3. Background field 

As the second major ingredient of a more realistic crust 
model we have to specify background field profiles more 
justified with a view to NS physics than the ad-hoc as- 
sumptions employed in our former papers. That's why 
we turn at this point to a model simulating the Ohmic 
decay of the magnetic field in a cooling NS's crust 
(i.e., in a spherical shell), which is just the situation 
where we expect the most significant effects of the Hall- 
instability to occur. When thinking about initial con- 
ditions for such simulations, one has to cope with the 
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Fig. 4. Radial profiles of diffusivity r\ (solid) for different 
ages of the NS, and Hall coefficient a (broken). Top: FP, 
bottom: PS model. All r\ profiles for ages > 10 7 yrs (FP) 
and > 10 6 yrs (PS) practically coincide. 



fact that, unfortunately, there is less certainty about the 
very origin of the NSs' magnetic fields. Favorite mecha- 
nisms of magnetic field gener ation are proto-NS dynamos 
({Thompson fc Duncanlll993t) and thermo electric instabil- 
ities (see, e.g.. IWiebicke fc Geppertlll99fil) . but inheriting 
the magnetic field from the NS's progenitor seems well 
possible, too. Since the way of field generation is surely 
significant with respect to the details of the field struc- 
ture, there seems to be no other resort when having to 
specify it than the principle of greatest possible simplic- 
ity. Therefore, we decided to consider dipolar fields only. 
Further on, it seems to be reasonable to assume with re- 
spect to their radial profiles that not all parts of the (in the 
case of a proto-NS dynamo: later) crust equally partici- 
pate in the field generation. Another process to be taken 
into account, is the fallback accretion burying the mag- 
netic field again within the crust soon after its emergence. 

On the other hand and somewhat conflicting with max- 
imum simplicity, we had to ensure that enough curvature 
is 'i njected' into the radial magnet ic field profile initially. 
See iRheinhardt fc Oepperd <j2002h for the essential role 
of the second derivative of the background field's profile 



(which is partly taken over by the gradient of a in the 
present case as discussed above). 

At this point it is necessary to discuss the relationship 
between the only z-dependent Bo profile used in the plane 
slab geometry, and a both r- and ^-depending profile re- 
sulting from decay calculations in a spherical shell. Since 
£? as assumed in (0 must be parallel to the boundaries, 
the plane slab model is only meaningful in the close vicin- 
ity of the magnetic equator 9 = %/2 (9 - polar angle). 
f(z) for a certain instant r has therefore to be identi- 
fied with the ^-component of B(r, 9, r) where r = R — z, 
= tt/2. Accordingly, the initial background field is com- 
pletely specified in geometry and strength after having 
prescribed f(z) for r = 0, which we denote as / in . (The 
restriction to background fields parallel to the slab bound- 
aries can be relaxed to fields with a constant vertical com- 
ponent. We will treat this case in a forthcoming paper.) 

In view of all these aspects, we generally chose initial 
profiles showing at least one maximum inside the crust 
and being zero beneath a certain initial penetration depth 
z- m . We fixed the latter by prescribing the corresponding 
density p- in to be either 10 13 or 10 12 gcm~ 3 . The smaller p m 
results in larger derivatives of the background field profiles 
thus possibly favoring the instability at early instances. 
On the other hand the decay will be accelerated and we 
have to expect smaller growth rates at higher ages (see 
Sect. I3J) . The following initial profiles of the 3rd and higher 
degrees were chosen: 



fin(z) 
fin(z) 



-2B n 



R 



R- z 



i?\ 4 z 3 



z- m J R 2 {R-z 



(1- 
■(1 



) cubic (13) 



) heptic (14) 



R R ( z 
nBd sin I 4tt — ) sinusoidal (15) 

Zin Z V Zi, 



Here, is just the initial polar surface magnetic field. 
Note, that the denominator R—z in (|13|l (|15[) doesn't in- 
fluence the degree of f- m with respect to z significantly 
since inside the crust R—z R. 

The decay calcula tions are straightforward (for details 
see IPage et al.l l200f]|) : we only mention that the back- 
reaction of the magnetic field on the temperature distribu- 
tion and cooling history via Joule heating was neglected. 
Hence, the results depend linearly on the initial condi- 
tions. The same boundary conditions as referred to for 
the magnetic field perturbations were used. However, we 
refrained from adjusting the initial profiles (|13|) - l|15|) to 
the vacuum boundary condition at the surface since the 
decay simulation will do this job automatically. (Likewise, 
diffusion smoothes out the 'kink' at z = z- m .) All runs were 
performed with B& = 10 13 G and snapshots of the decay- 
ing field were taken at nine instants from 3 x 10 3 to 3 x 10 7 
yrs in order to define the profiles f(z). Additionally, we 
scaled them (that is, B^) by factors 2 or 5 to study the 
dependence of the instability on the initial field strength. 

Figures El to show typical examples of the result- 
ing f(z) for both models. (See also Fig. ^ for schematic 
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sketches of characteristic background fields.) Moreover, 
these figures present the magnetization parameter wg r e 
and — only for the FP model, in Fig. [3] — the quantity 
( a f'Y / r h further on referred to as "curvature parameter". 



p (g/ccrri) 




1.0x10 2.0x10 3.0x10 4.0x10 
depth 2 (cm) 

p (g/ccm) 




1.0x10 2.0x10 J. 0x1 4.0x10 
depth 2 (cm) 



p (g/ccm) 
2.8x10 12 3.0x10' 




1.0x10 2.0x10 3.0x10 4.0x10 

depth 2 (cm) 

p (g/ccm) 

1.0X10 10 1.4X10 11 2.8X10 12 3.0X10 13 
i 1 1 1 

3.e + 03 yrs 
l.e+04 yrs 
3.e+04 yrs 
l.e+05 yrs 
3.e+05 yrs 
l.e + 06 yrs 




1.0x10 2.0x10 3.0x10^ 4.0x10 
depth 2 (cm) 

p (g/ccm) 

1.0x10 10 1.4X10 11 2.8x10 12 3.0x10 13 
8x1 0~ b 




1.0X10 2.0x10 3.0x10 4.0x10^ 
depth z (cm) 



Fig. 5. Background field Bo, magnetization parameter u)B a Tc 
and curvature parameter {aB' )' /rj resulting from the FP 
model with sinusoidal initial profile 1131 and initial penetra- 
tion density p\ n = 10 13 gcm -3 , and Bd = 10 13 G. Red parts 
indicate adherence to the sign condition 1161 . 
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Fig. 6. Background field Bq and magnetization parameter 
UB T e resulting from the PS model with cubic initial pro- 
file (|13J) and initial penetration density p- m = 10 13 gcm -3 , 
and B d = 10 13 G. 



It is an estimate for the ratio of the potentially energy- 
delivering term curl(a curl Bq x b) to the (anyway) energy- 
dissipating term curl (n curl b) in 10 . (Note that the term 
curl(acurl6 x Bq) is energetically neutral.) Therefore we 
suppose that the curvature parameter rather than u)B T e 
itself is the key indicator with respect to the vigour of the 
instability. For the PS model, the curvature parameter is 
only piecewise continuous due to the non-smooth profiles 
of 77 which result in /" profiles only piecewise continuous, 
too. Nevertheless, the solutions of Q are well defined, at 
least in the mathematically weak sense. 



Fig. 7. As Fig. but with sinusoidal initial profile (|15fl . 



On the other hand, we want to stress that the magni- 
tude of the curvature parameter alone is again not suffi- 
cient to infer the existence of the instability, which relies 
on a complicated interplay of vector quantities. Therefore, 
an additional qualitative criterion must be satisfied in or- 
der to have the potentially energy-delivering term in J3J 
really energy-delivering. From simple qualitative consider- 
ations with axisymmetric fields in a homogeneous medium 
it is possible to conclude that f"{z)f{z) < has to be 
satisfied, at least over a certain z-interval. Or generalized: 
Background field and curvature parameter must have dif- 
ferent signs: 

(af'y(z)f(z) < . (16) 
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3. Results 

3.1. Numerics 

The system of ordinary differential equations @ was dis- 
cretized by help of symmetric difference formulae of sec- 
ond order; near boundaries unsymmetric formulae were 
employed, if necessary. In most cases an equidistant grid 
with a typical node number of 1200 was used, but in quite 
a number of spot checks of convergence the node number 
was repeatedly doubled up to 9600 maximum. 

For the numerical solution of the resulting com- 
plex non-Hcrmitian standard matrix eigenvalue prob- 
lem we took benefit of the package ARPACK (routines 
znaupd , zneupd) . Some results were checked applying the 
LAPACK routine zgeev. The results from both packages 
agree in at least six digits (but comparisons were possible 
for moderate node numbers only). 

Unfortunately, the order of the difference formulae was 
not always reflected by the convergence rate of the eigen- 
values with respect to the node number: For the PS model 
convergence was usually significantly slower. Perhaps this 
behavior is connected with the discontinuities of the co- 
efficients rf , a' and /" described above. This supposition 
is supported by the fact that the convergence rates of the 
FP model are close to (although not exactly) quadratic. 

3.2. General aspects of the results 

As with the homogeneous density model, both oscillat- 
ing and non-oscillating unstable modes exist, where the 
fastest growing mode of any specific model turned out to 
be always amongst the latter. As another common feature 
of the fastest growing modes, their wavenumber k y has 
always been found to be zero. When deriving the back- 
ground field from an axisymmetric poloidal one in the 
vicinity of its magnetic equator (as we did), association 
of the y-co-ordinate with longitude is surely the proper 
choice. Thus, with some care one may suppose that in a 
spherical shell the most unstable modes are preferentially 
axisymmetric. 

Figure 00 shows growth rate 3?(p) and oscillation fre- 
quency 3(p) as functions of k x and k y for a typical ex- 
ample (PS model, age = 3 x 10 5 yrs, cubic initial field, 
Bd = 5 x 10 13 G , pin = 10 13 gcm -3 ). The rectangle enclos- 
ing the unstable region in the fc^-fcy-plane is defined by 
OA/d < k x < 18.4/d, < k y < 34.4/d But, when obey- 
ing the background dynamics permissibility condition (see 
Sect. 12 only a considerably smaller unstable region can 
be considered feasible. 

In contrast to the growth rates, the oscillation fre- 
quencies grow in general with k y ; their highest values 
seem to occur always at the edge of the unstable re- 
gion in the fc^-fcy-plane, that is, for 3?(p) — » (and 
k y ^ 0). Note, that the oscillating unstable modes can 
be r egarded as helicoidal wa ves with growing amplitudes. 
(Cf. IVainshtein et al.l i|2000l) who considered damped he- 
licoidal waves in a stratified crust.) 



In Figs. El to E| the growth times of the fastest grow- 
ing modes l/max(3?(p)) (simply referred to as "growth 
time" Tg r owth of a specific model) for the most important 
cases considered, are presented as functions of age and 
initial field strength Bd- In addition, the tangential pe- 
riod lengths of these modes, A max , are given; because of 
£,max _ q q Uan tity is defined as 



It should be regarded as one of the dominating scales of 
the unstable modes since, in contrast to the homogeneous 
density model, in the majority of cases no prominently 
small radial scales were found (see Fig. I15p . Instead, the 
radial scales of the unstable modes seem to be determined 
simply by the radial scale of the background field. Except 
for some earlier stages of the PS model, A max depends in 
general only slightly on Bd- Since the wavenumber k x was 
given only discrete equidistant values, the period lengths 
for different values of Bd frequently even coincide (indi- 
cated in the figures by filled symbols instead of open ones) . 

3.3. Growth times 

We stress again that the growth times have to be consid- 
ered as referring to snapshots only. That is, the value for 
a given NS age has been calculated assuming that the in- 
stability starts just at that age, with the background field 
given at exactly that moment. Clearly, the results must 
not be interpreted as a sort of temporal evolution of the 
growth times: Any occurrence of the Hall-instability will 
affect the strength and structure of the background field, 
thereby changing the conditions for the occurrence of the 
instability itself at later moments. This shortcoming of the 
results induces us to emphasize once more the necessity of 
full non-linear calculations. 

As an overall property, we state that the dependence of 
Tgrowth upon the initial polar magnetic field Bd is always 
monotonically falling. Hence, we never entered the range 
of B ( ] where higher values may y ield larger growth times 
(see iRheinhardt fc Gepperdl2002|) . 

3.3.1. FP model 

The growth times of the FP model (Figs. El to EJ Tab.EJ 

have in common that they start with small values (most 
in the order of magnitude 10 3 yrs and smaller) at the 
youngest age, exhibit a nearly linear dependence on age 
(a power law with an exponent between 0.8 and 1) until 
10 5 ... 3 x 10 5 yrs, reach a maximum at 10 6 . . . 3 x 10 6 yrs, 
and fall progressively with age later on. At early stages 
the background dynamics permissibility condition is well 
fulfilled for Bd = 5 x 10 13 G but must be put in ques- 
tion in most cases with Bd = 10 13 G. Analogously, the 
maximum growth times are reliable for most cases with 
Phi = 10 13 gcm~ 3 , but questionable for some of the cases 
with p in = 10 12 gcm~ 3 and B d = 10 13 G. At the latest 
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Fig. 8. Eigenvalue p as function of the wave numbers k x , k y for the PS model with cubic initial background field profile 
<pH5|) . initial field B<± = 5 x 10 13 G, initial penetration density p ln = 10 13 gcm -3 , age 3 x 10 5 yrs. Left: growth rate 
negative values were set to zero; right: frequency 3(p), not shown for 3?(p) < 0. p in units of 4/(<f/cm) 2 s _1 = 
2.68 x 10 _11 s _1 , wave numbers in units of 2/d = 5.18 x 10 _6 cni -1 . 



stage considered all the growth times satisfy the condi- 
tion. We note that the dependence of the growth time 
on £?d is not far from being linear for p ln = 10 13 gcm -3 , 
whereas for p m = 10 12 gcm -3 the dependence is of higher 
order. In this case, with the cubic profile and Bd = 10 13 G 
even a "gap" in the growth time curve occurs: From the 
age of 10 4 to the age of 3 x 10 5 years there were no growing 
modes to be found at all. 



' I '\ 1 000 




10 3 10 4 10 5 10 6 10 7 10 1 



Comparing the growth times for different values of 
Pi n , one has to state that the smaller initial penetration 
depth (pin = 10 12 gcm -3 ), although being connected with 
stronger gradients of the background field, is nevertheless 
disfavoring the instability. This tendency becomes increas- 
ingly apparent with growing age. We explain it with the 
accelerated decay of these "shallower" profiles. 



Comparing the growth times for different initial profile 
types, one can see that the values for the heptic and the 
sinusoidal profiles are close together at the earlier stages 



Fig. 9. Growth times and tangential period lengths of the 
fastest growing modes for the FP model with cubic ini- 
tial background field profile l|13|) . initial penetration den- 
sity p in = 10 13 gcm -3 . Solid, diamonds: Bd = 5x 10 13 G; 
dashed, triangles: Bd = 10 13 G. Thick lines, big symbols: 
growth times; thin lines, small symbols: period lengths. 
Full small diamonds: coincidence of the period lengths for 
different values of Bd- 
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Fig. 12. Growth times and tangential period lengths of 
the fastest growing modes for the PS model with cubic ini- 
tial background field profile l|13|). initial penetration den- 
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Fig. 10. As Fig.|Hl but with heptic initial background field 
profile lfH|) . 
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Fig. 11. As Fig. El but with sinusoidal initial background 
field profile fTK|) . 
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Fig. 13. As before, but with heptic initial background field 
profile |H|). 



and clearly smaller than the ones for the cubic profile. 
However, the late stages reveal the sinusoidal profile to be 
the favorite; in the case of p m = 10 13 gcm~ 3 it yields the 
smallest late stage growth time at all (20,240 yrs). As a 
rule of thumb, the growth times for the cubic, heptic and 
sinusoidal profiles behave like 4:2:1 at late stages. These 




10 3 10 4 10 5 10 6 10 7 10 8 
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Fig. 14. As before, but with sinusoidal initial background 
field profile (|T5|) . 



(G) 



profile type 


cubic 


heptic 


sinusoidal 


n global min. 


3,500 


2,920 


1,740 


3 global max. 


227,490 


71,230 


55,330 


x occurs at age 


10 s 


10 s 


3 x 10 7 


10 val. at max. age 


135,300 


69,630 


55,330 


global min. 


590,900 


68,590 


248,650 


occurs at age 


3 x 10 6 


3 x 10 3 


3 x 10 6 


S global max. 


3,071,190 


1,796,610 


4,540,120 


o 

1-1 occurs at age 


3 x 10 3 


10 4 


10 s 


val. at max. age 


658,290 


342,250 


271,080 


gap 


10 4 . . . 10 5 3 x 10 4 . . . 10 5 


10 4 . . . 3 x 10 4 


Table 3. Growth times r gr owth 


(yrs) of the PS model. 


For B d = 5x 10 13 


G the global 


minimum 


occurs always 



at the age of 3 x 10 3 yrs. 



ratios reflect even more or less quantitatively the relation- 
ship of the corresponding curvature parameters at the late 
stages, sec Fig. |SJ 

3.3.2. PS model 

A clearly pronounced and characteristic common property 
of the data shown in Figs.^Jto^J is the independence of 
Tgrowth on age for stages older than 3 x 10 5 yrs. Further on, 
for these ages the dependence of T gr0 wth on is linear to 
a high accuracy. We explain the first effect by referring to 
the weak dependence of the curvature parameter (ocf')'/r] 
on age (see Figs. |H]and[7Jl and the second by stating the 
dominance of the Hall terms in in this range of ages. 

In contrast, for the earlier stages the sensitivity of the 
growth times with respect to is extraordinarily high. 
It is highest for the cubic profile for which at the age of 
3 x 10 3 yrs a scaling of by 5 reduces the growth time by 
almost 3 orders of magnitude. In general, with falling a 
maximum of the growth time is forming between the ages 
of 10 4 and 10 5 yrs. "Gaps" in age without any unstable 
modes are listed in Table |3 
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B A (G) 




profile type 


cubic 


heptic 


sinusoidal 




p in (gcm" 3 ) 


10 13 


10 12 


10 13 


10 12 


10 13 


10 12 


<n 

rH 

o 

rH 

x 
in 


global minimum 
global maximum 
value at highest age 


OOU 

256,060 
84,880 


i 7nn 
1, lUU 

1,150,350 
326,520 


Kin 
OlU 

107,650 
38,780 


561,050 
185,630 


ooO 
53,590 
20,240 


con 
63U 

382,800 
122,370 


o 

rH 


global minimum 
global maximum 
occurs at age 
value at highest age 


7,880 
1,280,780 
3 x 10 6 
358,870 


74,900 
12,767,800 
10 6 

1,466,830 


3,470 
536,040 
3 x 10 6 
177,270 


8,290 
3,711,550 
3 x 10 5 
720,110 


2,870 
266,740 
3 x 10 6 
88,050 


5,620 
192,520 
3 x 10 6 
517,920 



Table 2. Growth times T growt h (yrs) of the FP model. For B^ = 5 x 10 13 G the global minimum and maximum occur 
always at the ages of 3 x 10 3 and 3 x 10 6 yrs, resp. For B^ = 10 13 G the global minimum occurs always at the age of 
3 x 10 3 yrs. For the cubic profile and p; n = 10 12 gcm~ 3 no unstable modes exist between the ages 10 4 and 3 x 10 5 yrs. 



FP PS 

dependence on age falling constant 
during rapid cooling- 
age at which r gr0 wth some 10 6 some 10 s 
is maximum (yrs) 

minimum r growth in » 20,000 55,000 
late stages (yrs) 

sensitivity of T gr0 wth linear to extremely nonlinear 

with respect to B& weakly nonlinear at early ages, 

linear from 10 yrs on 



An early onset of the instability could have a noticeable 
impact on the distribution of the magnetic field strengths 
across the NS population. Namely, if we had a growth time 
of, say, 500 yrs at the age of 3,000 yrs changing to, say, 
1,000 yrs at the age of 10,000 yrs one would expect an es- 
sential part of the background field energy to be consumed 
by the unstable modes during some 10,000 yrs after birth. 
XXX However, we see at least two reasons to suggest a 
very careful use of these "early" growth times. First, dur- 
ing the early stages the instability may benefit from per- 
haps unphysically high values of the curvature parameter 
inherited from our initial background fields. For this, the 
value of the curvature parameter reached even infinity at 
P = Pin because the background field profile f(z) shows 
a kink there. For the later stages, instead, we are opti- 
mistic that this heritage lost its significance due to the 
long lasting period of smoothing diffusion. Second, even if 
the value of the curvature parameter could be regarded as 
realistic, one had to take account of the high sensitivity of 
the conductivity with respect to temperature during the 
hot young stage of the NS. Likewise, the high sensitivity 
of 7"g rowt h with respect to B^ and, hence, to the curvature 
parameter depending on conductivity, should be consid- 
ered. Thus, due to the enhanced local heating by the "hot 
spots" of the unstable modes (see Fig.ll7|l an efficient self- 
limitation of the mode amplitudes or even a switching-off 
of the instability may take place. Perhaps, this happens far 
before a significant part of the background field's energy 
is transferred to unstable modes. 

3.4. Eigenfunctions 

Figure 1151 shows some typical eigenfunctions of the 
poloidal and toroidal scalars, s and t, for the most favor- 
able background field profiles, that is, the sinusoidal ones 
with pin = 10 13 gcm -3 , £>d = 5 x 10 13 G at the earliest and 
latest stages, respectively. Table lists major qualitative 
differences of the eigenfunctions for the two NS models. 



Table 4. Comparison of the results for T grow th of the FP 
and the PS models 



The background dynamics permissibility condition is 
not well satisfied at all the young stages, except for B& = 
5 x 10 13 G. As the background field and diffusivity profiles 
don't change significantly for ages larger than 10 6 yrs, the 
growth times are reliable from this age on. A comparison 
of the different profile types with respect to the later stages 
yields the sinusoidal one to be the most favorable for the 
instability, with the heptic one being close on its heels. 

3.3.3. Comparison FP — PS model 

From a final view on the results for both models we sum- 
marize the main differences between them as shown in 
Table H 

As a striking common feature of a majority of all the 
specific models one finds the very short growth times 
at the youngest ages. T his is in contradiction t o our 
assumption, introduced in lGeppert fc Rheinhardtl l)2002j) 
and based on the magnitude of lobt c , that the Hall insta- 
bility switches on only after a considerable cooling down 
of the NS, i.e., only after some 10 4 to some 10 5 yrs. XXX 
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Fig. 15. Moduli of the scalars s and t of the fastest growing mode for the FP (left) and PS (right) model with the 
sinusoidal initial background field profile 1151) , initial field B^ 
at the ages of 3 x 10 3 yrs (above) and 3 x 10 7 yrs (below). 



5 x 10 G, initial penetration density pi n = 10 gem 







FP model 


PS model 


primary 
of \s\ and 


maximum 

t\ 


almost fixed near z/d = 0.4 except at last age 

when closer to surface: 

at z/d = 0.2 for |s| and z/d = 0.35 for \t\ 


almost fixed near z/d = 0.2 


secondary 


maxima of t| 


lower than corresponding primary maximum 
except at last age; 


always lower than primary maximum 

"early" secondary maximum at z/d = 0.1 at 
age 3 x 10 3 yrs 


inner zero 


of \t\ 


emerges at age 3 x 10 6 yrs, moves from 
z/d — 0.2 to z/d = 0.05 with growing age 


emerges at age 3 x 10 yrs, moves from 

z/d — 0.07 to z/d = 0.04 with growing age; 

"early" zero at z/d = 0.14 at age 3 x 10 3 yrs. 



Table 5. Comparison of qualitative features of the eigenfunctions (cf. Figs. I15land ll6(l . 



In general all eigenfunctions show one pronounced 
(primary) maximum in both |s| and \t\ with two addi- 
tional (secondary) local maxima of \t\ close to the vacuum 
boundary for the two (FP) and five (PS) last ages, re- 
spectively. The global maxima of |s| lie a bit closer to the 
vacuum boundary than the corresponding ones of \t\. 

With one exception, the secondary maxima are lower 
than the primary ones. They are sharpening with growing 
age and B^, and emerge immediately below and above a 
zero of t which exists from the age of 3 x 10 5 yrs on and 



moves towards the vacuum boundary with aging. In order 
to discuss this feature, we stress at first that because of its 
clear convergent behavior with respect to grid size it has to 
be considered as a physical fact rather than a numerical ar- 
tifact. Since the zero of t mimics a vacuum boundary (see 
(|10Jl h it is intriguing to relate the peak beneath it to the 
similar peak right beneath the v acuum boundary in the 
homo geneous density model (see iRheinhardt fc Genperd 
12002ft . 



14 



Rheinhardt, Konenkov & Geppert: Hall-Instability in Neutron stars 



Obviously, the perturbation fields concentrate mainly 
in regions of high values of -Bo and/or the curvature pa- 
rameter, respectively. 

To judge whether the field structures of the eigcn- 
modcs might be considered small-scale, we state first that 
their major radial scale corresponds with the major ra- 
dial scale of the corresponding background field. The tan- 
gential scales can simply be identified with (half of the) 
tangential period lengths A, which again reflect the ra- 
dial background field scale. When comparing with the 
crust thickness, there are surely no small scales, except 
the structures around the secondary maxima. But, when 
comparing with the NS perimeter it is well justified to 
claim small-scaleness since for the FP model the A max lie 
between 700 and 1000 m, whereas the NS's perimeter is 
67 km. For the PS model we have A max = 1, 500 ... 4, 500 
m with a perimeter of 103 km. 

Because the non-smoothness of the conductivity pro- 
files in the PS model is not directly expressing itself in 
the background field profiles, one could expect the same 
for the eigenfunctions. Indeed, the eigenfunctions of the 
first four ages don't show such signatures, but those of 
the five last ages do. This fact is only partly explainable 
by the sharper kinks in the r\ profiles of the last four ages 
in comparison with those of the first five (see Fig. . It 
strengthens again the impression that there is a qualita- 
tive difference between 'earlier' and 'later' stages of the 
PS model. 

Figure 1161 shows perturbation field structures for se- 
lected cases. As a general tendency the fields of the FP 
model reach deeper layers of the crust and are less con- 
centrated in z-direction than those of the PS model. At 
early stages the fields tend to be relatively stronger out- 
side the slab in comparison with their maxima inside. Note 
the very sharp maxima of b y near the surface for the lat- 
est stages of the FP model, which are connected with the 
occurrence of secondary maxima in the toroidal scalar t 
(see Fig. [To)). 

4. Summary and Conclusions 

The results presented here prove that the instability, origi- 
nally found in a slab of uniform density, does persist when 
a typical stratification of density and chemical composi- 
tion as present in NS crusts, and a realistic background 
field are employed. 

We found that independent of the specific NS-model, 
the strength and structure of the background field, the 
growth times are smallest at early stages of the NS, reach- 
ing then a maximum at a model depending intermediate 
age and become smaller again in the process of further 
cooling. In some cases (e.g., PS model, £>d = 10 13 G) 
the instability disappears for too low background field 
strengths in a period of intermediate ages. 

In order to understand this behavior, note that during 
the early period the crustal crystal is so hot that a large 
number of phonons is excited. Hence, the electron relax- 
ation time r e is relatively small, causing in turn a small 



magnetization parameter lobt e . But, competing with that, 
the magnetic field, initially confined to a shallow crustal 
layer had not yet time enough neither to decay nor to 
diffuse deeper into the crust therefore still preserving the 
curvature of its initial profile. It seems, that this effect 
dominates the relative smallness of u>BT e - In turn, during 
late, cool stages r e is much larger but the curvature of the 
background field profile is due to its progressing diffusion 
and decay considerably smaller. (See Figs. to 0). 

Quite general, ou r conclusion in 

iRheinhardt fe GeppertJ l|2002|) that the appearance 
of the instability depends both on a sufficiently strong 
curvature of the background field profile and a sufficient 
large magnetization parameter wsr e has thus been 
confirmed. However, it turned out that in models with 
stratification one has to include the derivative of the 
Hall-coefficient into the notion of 'curvature' by introduc- 
ing a properly defined curvature parameter (see Sect. 2). 
It estimates the ratio of the most relevant terms in the 
governing equations © ■ 

The differences between growth times within the set of 
early stages (instants labelled 1 to 5 in Figs. |51to[7|) or 
late stages (instants labelled 7 to 9) , respectively, are well 
reflected by the curvature parameter. But, unfortunately, 
it fails when comparing growth times between the two sets 
as it predicts the shortest growth rates for the late stages. 
In both NS models studied, there seem to exist subtle and 
hidden features of the coefficients and background field 
profiles which cause a qualitative contrast between early 
and late stages not understood up to now. 

Our results give rise to the idea that the Hall- 
instability will act relatively early in the NS's life. Hence, 
it could reduce and smooth out the background field in its 
progress, so that later on (say after 10 5 yrs) the conditions 
for the occurrence of the instability are no longer given. 

On the other hand we found (mainly for the PS model) 
a high sensitivity of the 'early' growth times with respect 
to the background field strength. An increase by a fac- 
tor of 5 can cause a decrease of the growth time by al- 
most three orders of magnitude. Recalling the importance 
of the profile curvature, we conclude that details of the 
background field (strength and structure) determined by 
the generation process of the NS magnetic field may play 
a crucial role with respect to appearance and vigour of 
the instability at early stages. Since these processes are 
far from being completely understood up to now, the only 
detail one can discuss is the initial radial extent of the 
generated fields. Generation by a thermoelectric instabil- 
ity in the relatively thin liquid shell forming the later 
crust surely provides radial profiles with a lot of curva- 
ture. But whether they arc really instability-friendly has 
to be put in question as a small initial penetration depth 
of the background field (p ln = 10 12 gcm -3 ) turned out to 
be non-favorable. Considering alternatively ty pical struc- 
tures of convection in a n ascent NS (see, e.g.. iKeil et alJ 
Il996t iFrver fc Hegel l2000|) one might conclude, that the 
assumption of bigger initial penetration depths could be 
appropriate for magnetic fields generated by a proto-NS 
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dynamo. In all, we suggest to make use of the 'early' 
growth times with great care only. 

The late stages, in contrast, are increasingly less influ- 
enced by peculiarities of the somewhat arbitrary assump- 
tions on the initial fields. At these stages the dependence 
of the growth times on the background field strength is 
undramatic (close to linear). So we think that we can rely 
safely on the growth times obtained at ages > 10 5 yrs. 

If the Hall-instability sets on at late stages, small- 
scale perturbations amplified at the expense of the back- 
ground field may survive for a long time. To estimate 
it, let us first look onto the Ohmic decay times for typi- 
cal perturbations. Their scale lengths are in the order of 
5 x 10 4 . . . 10 5 cm, and their field maintaining currents cir- 
culate in a depth of 2 ... 5 x 10 4 cm corresponding to a 
density of about 10 12 g cm -3 . Taking into account that 
at an age of 10 5 . . . 10 6 years the conductivity in that re- 
gion is in the order of 10 26 . . . 10 27 s" 1 , we arrive at values 
of 10 7 ...10 8 years. XXX However, these perturbations to- 
gether with the background field are subject to the Hall 
effect. Following a usual argumentation, one would there- 
fore expect that the real decay times are shorter than the 
ohmic ones. This conclusion would be based on the idea, 
that the Hall effect, although being conservative on its 
own, leads to an accelerated decay of a large scale field 
because energy is flowing along the Hall cascade towards 
small scales the dissipation rate of which is higher than 
that of the large scale field. However, this is not the whole 
truth for a situation in which energy is already concen- 
trated in a small scale. Along with a Hall cascade to 
even smaller scales dissipating more quickly another one 
to larger scales can occur which dissipate more slowly. 
Whether or not such a double-sided cascade leads to an 
overall acceleration of the decay of the original small-scale 
mode must remain open. XXX 

The possible occurrence of relatively persistent small- 
scale field structures at the NS surface together with the 
(episodically) accelerated decay of the large-scale back- 
ground field caused by the energy transfer to small-scale 
modes during their rapid growth represent the most im- 
portant aspects of the Hall-instability with respect to ob- 
servations. Possible observable consequences arc therefore 

— deceleration of the NS spin-down due to the accelerated 
magnetic field decay. XXX Note, that the real extent 
and duration of the latter can only be determined by 
analyzing the nonlinear stage of the instability includ- 
ing saturation. XXX 

— a hotter NS surface as a consequence of enhanced Joule 
heating due to the concentration of magnetic energy in 
the small scales of unstable eigenmodes 

— glitches and bursts as a consequence of enhanced 
Lorentz forces making crust cracks more probable 

— emergence of radio subpulses due to the existence of 
small-scale structures in the vicinity of the magnetic 
pole. 

A discussion of the first thr ee effects can be found in 
iGepnert fc Rheinhardtl 1 20021) . We will not repeat it here 




Fig. 17. Perturbation of the Joule heat sources density 
(77/271") curl So • curl 6 in arbitrary units corresponding to 
the perturbation field of the left column, third row of 
Fig. E| Green to red — positive, green to blue — negative 
deviations from the background heat sources. 

since all the general arguments hold true for the Hall- 
instability in a stratified crust, too. However, some modi- 
fication seems to be necessary with respect to the enhance- 
ment of Joule heating, since in our present results the case 
of "hot spots" close to the surface is not the prevailing one 
(Fig. 117(1 . Thus, the argument, that small-scale tempera- 
ture features due to localized Joule heating will be washed 
out and remain unobservable in the light curve seems to 
be the more valid. On the other hand, an observable rise 
in the average surface temperature in comparison with 
the standard cooling due to the additional "deep" Joule 
heating is likely at least if the instability acts during the 
phase for which the standard model predicts a rapid cool- 
ing (10 5 . . . 10 6 yrs). 

A necessary ingredient of the pu lsar vacuum gap model 
of iRuderman & Sutherland! is the existence of rel- 

atively strong and small-scaled poloidal magnetic field 
structures just above the surface of the polar cap. In a 
forthcoming paper we will demonst rate that the cha rac- 
teristics of such a field as derived bv lGil et all j2003h can 
be found with some of the unstable modes presented here. 
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Fig. 16. Fields of the fastest growing modes for the FP (left) and PS (right) model, initial field = 5 x 10 13 G, initial 
penetration density p m = 10 13 gcm~ 3 . Upper two rows: cubic lfH5|) . lower two rows: sinusoidal lfT5|l initial background 
field profile. Ages as indicated in the middle (in yrs). The field component perpendicular to the paper plane is indicated 
via color encoding. 



